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ABSTRACT 

Shipboard measurements of the vertical profiles of mean 
humidity, temperature and wind speeds over the ocean have 
been investigated to verify their suitability to accepted 
turbulence similarity theory, In addition, the profiles have 
been used to obtain estimates of the temperature structure 
parameter. One hundred eighty-two profiles observed over 
a seven month period in 1974 aboard the R/V Acania in 
Monterey Bay have been compared to previous experimental 
results obtained in defining the meteorological parameters 
of friction velocity, roughness length, drag coefficient 
a.nd the Richardson number. The temperature structure param- 
eters computed from the profiles were provided to Lund (1975) 
for comparison with direct measurements. 

In general, favorable results were obtained in the param- 
eter comparisons which provided additional support to the 
suitability of shipboard profile measurements defining the 
marine surface layer. However, some evidence exists in 
roughness lengths and drag coefficient values to indicate 
that the ship's structure may be a factor influencing results 
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I. 



INTRODUCTION 



The air-sea interface region has many physical properties 
yet to be fully understood. Only in the past decade has 
there been a concerted effort to make measurements in the 
planetary boundary layer over the ocean and to examine exist- 
ing laws describing it. The Navy is actively engaged in 
marine atmospheric boundary layer research for many reasons. 

One of which is to accurately access the environmental 
influence on optical wave propagation over the ocean. Exper- 
imental evidence suggests that certain near surface atmospheric 
parameters can give reasonable results in specifying the 
optical conditions for the first tens of meters above the sea. 
In this area, problems involving applied theory, shipboard 
instrumentation and data processing exist. 

A primary purpose of this work is to evaluate, for the 
marine boundary layer, expressions which have been formulated 
on the basis of previous experiments. 

Data examined in the study were measured from the Naval 
Postgraduate School research ship, the R/V Acania shown in 
figure 1. Beginning with accepted basic surface layer theory 
and advancing to turbulence similarity theory, observed data 
is analyzed in part for its practicability in defining optical 
wave properties in an ocean environment. 

The suitability of the R/V Acania as a platform for ocean 
micrometeorological measurements was previously examined in 
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Research Vessel Acania 



Figure 1 
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detail by Cavanaugh (1974) and Welsh (1974) and verification 
of their work is a secondary objective. 

Examined in this study are ten minute time averaged 
values of wind, temperature and relative humidity measured 
at generally four levels. The observations were made during 
the months of February, March, August and September 1974 
while the R/V Acania was at anchor approximately two miles 
offshore in Monterey Bay, Coincident with these measurements 
were measurements to obtain spectra of the wind and temper- 
ature f lunctuations and also measurements of optical propaga- 
tion characteristics. The agreement of the profile data with 
boundary layer theory predictions and results of others will 
be examined to determine the validity of the shipboard 
measurements , 

The boundary layer parameters examined are the roughness 
length (z^). friction velocity (u^), drag coefficient (C^), 
Richardson number (Ri), and the temperature structure param- 
eter (C^). Examinations of some of these parameters are 
accomplished by applying the data to different calculation 
methods. 

The structure parameter is calculated from an indirect 
method based on similarity turbulence theory as proposed 
by Wyngaard, Izumi, and Collins (1971). The results were 
used in comparisons by Lund (1975), who analyzed direct 
measurements of this parameter for the same data periods. 
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II. THEORETICAL BACKGROUND 



A . GENERAL 

Investigations of temperature and wind fluctuations in 
the surface layer provide a basis for defining the stratum 
dynamics by turbulence similarity theory. Measurements of 
mean wind (u), mean temperature (T) and mean relative 
humidity (q) over a time scale greater than that of the 
fluctuating motion are used to estimate the vertical fluxes 
of momentum and heat due to the fluctuations, This can be 
accomplished because profiles of the aforementioned parameters 
are maintained by the momentum and heat fluxes. Turbulence 
similarity theory which is often referred to as the Monin- 
Obukhov theory relates the profiles to the fluctuation 
properties. Statistical analysis is used to arrive at 
empirical expressions which are accepted and are consistent. 
Dimensional analysis was used in the formulation of all 
expressions . 

The empirical relationships relating gradients of wind, 
temperature, and humidity to the two fluxes are 
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starred quantities of wind (u^) and temperature (T^) 
include the vertical flux parameters as indicated in the 
following expressions. 



. — 
= (-UW) 




(4) 

(5) 

( 6 ) 



B. PROFILES IDENTIFYING SURFACE LAYER CHARACTERISTICS 

Logarithmic plots of mean wind profiles can be used to 

define certain surface layer properties such as roughness 

length (z ), friction velocity (u*) and drag coefficient 
o ^ 

(Cjj). The latter is normally calculated for a ten meter 
height . 

Integrating equation (1) and assuming neutral or near 
neutral stability conditions, where 4)^(z/L) = 1, yields the 
equation 

z 

k z • (7) 

The wind speed at the surface, u(z ) is assumed to be zero 
for reasons of profile analysis. However, that may not be 
the case if the surface is also in motion, as is the case 
over the open ocean. Wind measurements taken at several 
levels and plotted as a logarithmic z profile yield an inter- 
cept of z^ and a slope of k/u,,,. Figure 2 is a representation 
of theoretical versus actual wind profiles in the marine 
surface boundary layer. 



13 



W- 






U Wind Speed 1m/sec 



Figure 2 



The drag coefficient is defined as the ratio of the squared 
friction velocity to the squared mean wind speed at a specific 
height . 
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C. DESCRIPTIONS OF STABILITY INFLUENCE ON THE SURFACE LAYER 

The functions (p , 4>„, are used to account for the 
^m' H’ ^E 

stability influence in equations (1), (2), and (3), The 
stability of the surface layer regime plays an important 
role in determining flux intensities, Monin-Obukhov theory 
determines the level where mechanical and thermal production 
of kinetic energy are equal to a dimensionable length (L) , 
The ratio of the height of the measurement (z) to the Monin- 
Obukhov length (L) serves as a stability index, z/L. The 
Monin-Obukhov length is defined as 



L = 



T u*' 

V * 



(9) 



where 

+ .61 q^ T (10) 

In unstable conditions the stability parameter (z/L) is 
approximately equal to the Richardson number (Ri) which is 
often used as the stability parameter, The following relation- 
ships between the stability parameter and the Richardson 
number have been proposed by Dyer and Hicks (1970) (unstable) 
and Webb (1970) (stable). 
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a is an empirically derived constant equal to ,5, 



D. INDIRECT CALCULATION OF THE TEMPERATURE STRUCTURE 
PARAMETER 



The temperature structure parameter (C^) is the most 



important variable in determining the optical propagation 
characteristics of the atmosphere, One means for calculating 
C^ is by employing turbulence similarity theory. Since 
atmospheric refractive index fluctuations are not included 
in these calcula.tions , the method indicated is an indirect 
estimation of atmospheric optical properties. In similarity 
theory estimates of the dissipation of turbulent kinetic 
energy (e) and dissipation of temperature variance (Gq) 
provide the link between direct methods such as proposed by 
Corrsin (1951), equation 13, and the conveniently relatable 
vertical fluctuations of wind and temperature in the surface 
layer . 



Research by Wyngaard, Izumi and Collins (1971) has resulted 
in nondimensionalized expressions relating atmospheric 
stability to the structure parameter. Dimensional analysis 
of the height (z), the friction velocity (u,^,), the mean 
potential temperature gradient produce 



4,^(k ) = (5 cej 
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(13) 
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functional expressions dependent only on the Richardson 
number. Similarity principals are elucidated by using the 
nondimensional stability parameter Richardson number defined 
as equation 14 for indirectly obtaining a comparable structure 
parameter to identify the optic properties existing in the 
surface layer. 



Ri = 









(14) 



Equations including the Richardson number functions 
(fi, fz, fa) are 
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The latter expression relates the temperature structure 
parameter and stability parameter and allows for an empirical 
determination of the fa function at varying conditions of 
stability. If the function fa is defined, the temperature 
structure parameter can be readily calculated from more easily 
measured parameters, and Ri, 
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III. EXPERIMENTAL OBSERVATIONS 



A. DATA COLLECTION 

The Naval Postgraduate School R/V Acania was the obser- 
vational platform. The ship was anchored in Monterey Bay 
at a location approximated in figure 3, This location 
facilitated complimentary and simultaneous optical experi- 
ments from two shore stations, Hopkins Marine Station and 
the Holiday Inn. Normally, observational periods were three 
days in duration beginning with a one day equipment calibra- 
tion and set up, continuing with a one day on station systems 
test and finishing with one day of data collection, Data 
collection was mainly during evening hours. The entire 
operation required considerable inter-departmental coordination 
between the meteorology and physics departments. 

Data was collected and analyzed for the dates and time 
shown in table 1. These observations represent over ninety 
percent of all the mean profiles obtained in Monterey Bay 
during the period from 25 February to 20 September 1974, 

Figure 4 shows the general instrument configuration 
used during the experiments. Variations in the arrangement 
occurred primarily because of individual component failures, 

B. EQUIPMENT 

The mean wind measurements were made with a Thornthwaite 
Associates cup anemometer wind profile register system. 
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MONTEREY scale, miles 



DATE 


TIME 


NUMBER 

OF 

LEVELS 


NUMBER 

OF 

PROFILES 


26 February 


1113-2225 


4 


23 


27 February 


1643-2150 


4 


25 


14 March 


1703-2138 


2 


20 


27 March 


1823-2142 


3 


17 


28 March 


1824-2150 


4 


14 


13 August 


1630-2200 


4 


26 


18 September 


1617-2202 


4 


24 


19 September 


1724-2353 


4 


25 


20 September 


0004-0644 


4 


8 



Table I 




Figure 4 
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model number 104. In operation the shaft of a three cup 
anemometer unit serves as the shutter between a light source 
and a photocell for each revolution. The cups are plastic 
cones reinforced with aluminum frames. They are attached 
to the rotating shaft by stainless steel tubing spaced at 
120 degree intervals about the shaft as shown in figure 5. 

The three cup assembly sets along with the other sensors 
were positioned at four levels on the bow with electrical 
leads to the after deck house laboratory. The sets have the 
characteristics of low starting ■ speeds with a small amount 
of internal friction which aids in checking inertial over- 
shoot. 

Temperature sensitive quartz crystal probes, (Hewlett 
Packard model HP-2850) were used to measure mean temperatures 
at each level, RF signals from the crystal probes and from 
a reference oscillator were mixed in the HP-2801A readout 
unit to produce a beat frequency whose signature can be 
analyzed to within 0.001 degrees centigrade per hertz. 

Each sensor simultaneously received pre-experiment calibra- 
tion against a platinum resistence wire thermometer in a 
temperature controlled circulating v/ater bath over the 
expected temperature range. The accuracy in achieving a 
.005 degree centigrade correction factor was a constant for 
each probe. A 3.7 meter flexible coaxial cable is permanently 
attached to the sensor head and the mast mounted probes are 
housed in an aspirated shelter as depicted in figure 6. 
Temperature values were automatically recorded on a printer 
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C,W. Thornthwaite Anemometer Cups 



Figure 5 
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Quartz Thermometer Probe 



Figure 6 
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tape. The Hygrodynamics Digital I hygrometer indicator and 
the Dunmore-type lithium chloride sensor were used to measure 
the relative humidity. The Dunmore-type lithium chloride 
sensor is illustrated in figure 6, The equipment operates 
on the basis of resistence change within an electrolytic 
solution generating a reference voltage variance which is 
proportional to the relative humidity change. Automatic 
temperature compensation in the instruments meet the follow- 
ing specifications for relative humidity, 

+ 3% relative humidity below 90% 

+ 4% relative humidity above 90% 

Calibration of the sensors was accomplished by comparative 
methods using saturated salt solutions in an enclosed 
chamber . 

Relative humidity was observed and recorded similar to 
that of temperature, as a printed output. A more detailed 
description of the devices, their calibration and deploy- 
ment can be obtained from Cavanaugh (1974), 
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IV, ANALYTICAL DATA PROCESSING 



A. DATA EDITING 

During the experiment, mean wind measurements over ten 
minute intervals were recorded in a log book. The recordings 
were the number of cup revolutions for the interval as 
recorded by electronic counters. Temperature readings were 
printed on a paper tape from a monitoring digital readout 
at approximately one second intervals, Initially relative 
humidity measurements were recorded the same as the tem- 
peratures. However, due to a printer malfunction, these 
readings were recorded by hand in log book after the 28 
March cruise. All information was later coded and key 
punched on computer cards for further analyses, 

Preliminary screening processes were performed on 
computer listings. The data were examined for gross errors 
or inconsistencies. Very few deletions were necessary at 
this phase. Next, individual mean wind and mean temperature 
measurements were plotted to give logarithmic graphical 
profiles for each time period. The computer program devised 
for this operation also performed a least squares fit of the 
data. 

In addition, mean wind profiles for certain periods 
were closely examined for any biases which might be caused 
by the ships' hull or superstructure. Approximately ten per 
cent of all the measurements was discarded either because 
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of incongruous profiles or uncomplementary data times. 
Rejected profiles were mainly those of wind.. They occurred 
principally when winds were very light. 

Temperature profiles were rejected when one or more of 
the automatically recorded readings appeared to be spurious. 
This was not discovered during the initial data inspection. 

B. DATA RECORDING ADJUSTMENTS 

Since the data was obtained over a seven month period, 
significant alterations in data recording took place due to 
equipment failure or experiment deviations. Mean wind read- 
ings were obtained at four levels for all the periods except 
when the number of levels were two and three for March 14 
and 27 respectively. In August, the order in which the mean 
wind levels were labeled changed from previous cruises. This 
data variation was compensated for by analyzing the measure- 
ments on a cruise to cruise basis with simple computer 
program modifications. There were no major problems encoun- 
tered in programmed analysis of the temperature data. When 
the indicated relative humidity was missing, adjustments were 
made by taking an average between the two adjacent available 
measurements . 

C. PROGRAMMING CONSIDERATIONS 

Cup anemometer readings were recorded in revolutions per 
time period. The relationship between revolutions to speeds 
in nonlinear. It changes with wind speed, A computer 
subroutine was devised to convert the number of cup rotations 
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into wind speed. The conversion accuracy for wind speeds 
decreases with decreasing revolutions, 

An IBM System/360 source library program called LSQPL2 , 
Least Squares Polynomial Fitting, was an integral part of the 
data analyses. LSQPS2 not only served to give linear func- 
tions for data profiles, but also was used as a means for 
approximating roughness length (z^) values, Once the least 
squares linear profiles had been calculated, profile values 
were estimated for levels one meter either side of the 
observation levels. This was done for the purpose of defin- 
ing gradients used in computations of the Richardson numbers. 

The basis for describing stability in discerning optical 
propagation conditions is the Richardson number. In testing 
for the feasibility of applying identical relationships over 
the ocean, Richardson numbers were determined from the 
Monterey data utilizing procedures similar to those set 
forth by Wyngaard et al . (1971), Mean profiles of temperature 
and wind were obtained by an analytical least squares fit of 
ten minute mean observations at their respective instrument 
heights. This provided the wind and potential temperature 
gradients necessary to compute a Richardson number at each of 
the four instrument levels. Equation 18 is the approximation 
for the mean wind speed gradient between two levels used in 
equation 14. 
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( 18 ) 



= ^ = u(z? ) - u(z, ) 
Az Z3 in 22/^1 



Z3 represents the geometric mean height, (Zj • z^)^, A 
similar expression was used to compute temperature gradient, 

^ . 

3 z 
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V. DATA ANALYZATION AND COMPARISONS 



A. FRICTION VELOCITY AND ROUGHNESS LENGTH RESULTS 

Important parameters describing the marine boundary 
layer are the friction velocity (u^) and roughness length 
(z^). The friction velocity is a measure of the momentum 
transfer described by equation 4 and the roughness length is 
a boundary parameter defined as the height where the loga- 
rithmic profile becomes a valid specification. These param- 
eters were calculated from the data by tv/o methods. One 
method utilized all profile levels, the other utilized only 
the wind speed at one level and the air water temperature 
difference. This was done to determine what difference would 
occur if the surface layer properties were obtained using 
data from one level instead of from a four level profile. 

The first method was based on a least squares fit of a 
stability corrected four level profile. The stability 
influence was determined at each level from the Richardson 
number. The following Richardson number categories were 
used in the corrections. 

-0,05 < Ri < ,001 neutral 

Ri > .001 stable 

Ri < -.05 unstable 

Stability relationships in equation 11 and 12 were used to 
obtain the stability length (z/L) which was then used to 
obtain the correction to the empirical expression for mean wind. 
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The plot of the mean wind versus the logarithm of 
height with the stability correction should theoretically 
give a straight line with a slope of k/u^ and an intercept 
of z^. Equation 19 was used to incorporate the stability 
influence into the profile 

in z - tp (^) = — u(z) + In z (19) 

The function ^(j-) was defined as follows from Paulson (1970) 

J-i 



iP(g) = 2 
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] - 2tan~^ x + ^ 
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- 16(g))^ 


unstable 


(21) 


^(f) 
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-4 7^ 


stable 


(22) 



A least squares fit was applied to the corrected profile to 
obtain the best estimate of the slope and intercept thus 
yielding the profile values for friction velocity and rough- 
ness length. 

The second method was to use a calculation suggested by 
Cardone (1969). He formulated an iterative procedure to 
find friction velocity and roughness length given a mean 
wind, an air-sea temperature difference and the instrument 
heights. By utilizing wind and temperature gradient expres- 
sions combined with equations similar to 9, an iterative 
algorithm was developed, which gives stability length 
convergence to a desired accuracy. The friction velocity 
was obtained using the stability criterion. The roughness 
length was then calculated from the friction velocity and 
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the following expression from Cardone, 



z = — + C,u. - C, 

O Z * 3 



( 23 ) 



The constants contained in equation 23 were derived so that 

the ten meter drag coefficient (C^q) would be a minimum 

when the mean wind speed was approximately six meters per 

second under neutral conditions. 

The results from the two methods are compared in figure 

7 which is a logarithmic plot of versus u^. Comparison 

is also made with Charnock's (1955) expression, 

2 



u. 



z = 



o a g a = 81,1 (24) 

The similar curvature of all three plots indicates agree- 
ment on the shape of the friction velocity and roughness 
length relationship. However, the least squares method 
produced roughness lengths on the order of one magnitude 
larger than Cardone 's values, 

In general, both methods appear to give good results, 
however, the average values obtained by the least squares 
method show less dispersion at low friction velocities than 
does Cardone 's one level method, 

A comparison of the roughness length versus the ten 
meter mean wind speed was made between results obtained by 
Ruggles (1970) and the least squares method as shown in 
figure 8, Ruggles hypothesised that the z^ peaks were 
associated with the effects of wind wave coupling. He 
designated the first peak as being coincident with the 
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FRICTION VELOCITY u^(cm/sec) 



Figure 7 
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Figure 8 
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minimum wind necessary for wave generation, The results of 
this study show similar maxirnums and minimums, most notibly 
the maximum near 3.5 meters and minimum near 4.5 meters, 



B. DRAG COEFFICIENT RESULTS 

Ten meter drag coefficients (C^^q) were calculated to 
provide further verification of the compatibility of the 
experimental profile results to results of others, Friction 
velocities which were computed by the least squares method 
and also using Cardone's expression and the ten meter mean 
wind (u^q) were used to obtain C^q by the expression 



u. 
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u 
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Plots of C^i^ versus u^ and C^^q 



(25) 

versus Richardson number 



appear in figures 9 and 10, Superimposed on the graphs are 
the results of others. 

In general, some aspects of the results compare favor- 
ably with results of others and indicate that over ocean 
shipboard application of turbulence similarity theory might 
be valid and feasible in the marine surface boundary layer. 

The least squares method produced average drag coefficients 
that were higher than the values of other experiments. This 
is attributed to the high values of friction velocity possibly 
due to the ship’s influence, The problem was addressed by 
Cavanaugh (1974), however, his results were not as pronounced 
as those in this study. 
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Figure 9 
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Figure 10 
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C. TEMPERATURE STRUCTURE PARAMETER 

The most important element for parameterizing optical 
propagation conditions in the marine surface boundary layer 
is the temperature structure parameter (C^), The indirect 
method was adopted to define at the four instrument levels 
where a computed Richardson number was available, First, it 
was necessary to use the Richardson numbers to define the 
fj function in equation 17 that had been empirically derived 
from the Kansas field data. Equation 17 was then used to 
obtain C^, Subsequently, the temperature structure parameter 
was transformed into the refractive index structure parameter, 
for additional comparisons with theoretical predictions, 
spectral and optical results, 

Calculated C^ decrease with height as expected for all 
stability conditions. Above a Richardson number of .20, the 
fj function in equation 17 is undefined. 

Temperature structure parameter (C,^) results were 
compared to direct measurement values by Lund (1975), 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The ship's motion and structure appeared to have had an 
effect on the wind profiles. This conclusion is based on the 
fact that the values of and C^q obtained by a least squares 
fit profile data were consistently higher than values observed 
in other experiments. 

Cavanaugh (1974) observed similar results and alluded 
to ship's motion influence on a logarithmic profile. Hull 
interference could disrupt the wind stream passing the 
instruments. This obstruction could cause the lower wind 
cups to slow and thus increase the apparent gradient. This 
would result in higher z^ values. This possible influence 
was examined in the analysis by disregarding the lower sensor 
level in the regression calculations. Roughness lengths were 
lowered slightly in most profiles examined, however, not 
enough to justify dropping this level in profile analyses. 

A suggestion for future experiments would be to position 
a sonic anemometer vertically over the side of the hull at 
main deck level to measure the forced vertical air motion 
caused by hull wind blockage. From data obtained in this 
manner, a ships wind stream signature for varying wind speeds 
might be found. This could find application in the more 
complex study of measuring turbulent motion while underway. 

With regard to temperature structure parameter estimates, 
the function f^ is required. Sufficient data is now becoming 
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available to permit evaluation of equation 17 for over the 
ocean observations. Comparison of the spectral results 
obtained by Lund (1975) and those computed in this study 
could accomplish part of this task. 

Because of the quantity of data to be screened only a 
cursory analysis of the observations has been completed in 
this study. More detailed examinations into the measurements 
would certainly yield additional knowledge, particularly with 
regard to the ship's motion and structure influences and 
instrument biases. 

It would be desirable to have mean wind observations 
closer to the ocean surface. Combined with the present 
profiles, a more detailed description of wind-wave coupling 
and its effect on turbulent flow could be undertaken, The 
additional measurements could be made from a buoy positioned 
upwind of the ship. Another- recommendation which is unrelated 
to application of turbulence similarity theory is that a small 
portable computer be obtained for future experiments to 
facilitate real time data analysis. The system now employed 
is extremely inefficient because many data discrepancies are 
not recognized until after the cruise. 
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